Engineering Quantum Jump Superoperators for Single Photon Detectors 
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We study the back-action of a single photon detector on the electromagnetic field upon a pho- 
todetection by considering a microscopic model in which the detector is constituted of a sensor and 
an amplification mechanism. Using the quantum trajectories approach we determine the Quantum 
Jump Superoperator (QJS) that describes the action of the detector on the field state immediately 
after the photocount. The resulting QJS consists of two parts: the bright counts term, represent- 
ing the real photoabsorptions, and the dark counts term, representing the amplification of intrinsic 
excitations inside the detector. First we compare our results for the counting rates to experimental 
data, showing a good agreement. Then we point out that by modifying the field frequency one can 
engineer the form of QJS, obtaining the QJS's proposed previously in an ad hoc manner. 

PACS numbers: 03.65.Ta, 42.50.Lc, 03.65.Yz, 42.50.Ct, 85.60.Gz 



I. INTRODUCTION 

Single photon detectors (SPD) represent the ultimate 
sensitivity limit for quantum photodetectors, and many 
quantum optics and quantum information applications 
are based on its existence Nowadays, there are sev- 
eral available types of SPD's sensitive to different light 
wavelengths and with a varying range of quantum ef- 
ficiencies @, S !, i, I, 0, S, 1 EE EL E2- Among 
many applications, SPD's are the main ingredient in the 
situations where one measures the electromagnetic field 
(EM) with few photons enclosed in a cavity, the pho- 
tons being counted one by one. Theoretical treatment 
of such a scheme, known as continuous photodetection 
model (CPM), was proposed by Srinivas and Davies (SD) 
in 1981 II 31 and did find many applications since then 

Hill El 

In SD scheme, in each infinitesimal time interval the 
photodetector has only two possible outcomes: either a 
single photon is detected ('click' of the detector), or it is 
not. In both cases the state of the field changes as time 
goes on: for a click, the field loses one photon and suffers 
a Quantum Jump; for no-click the field state modifies 
continuously and non-unitarily due to the monitoring of 
the detector [l?], Ell ■ Thus, besides allowing determina- 
tion of the statistical properties of the EM field through 
photocounts statistics, the detector also exerts a back- 
action on the field by which the outcome of the mea- 
surement modifies the cavity field state [l!|. This phe- 
nomenon was widely used for different theoretical pro- 
posals, e.g. for changing the field statistics from sub- 
Poissonian to super-Poissonian [2(J, for controlling the 
entanglement between two field modes (l5j or inducing 
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spin squeezing in a cavity [1 €31 ] - However, no experimen- 
tal verification of the CPM has been made until now, 
though nowadays it is quite realistic to make simple pho- 
tocounting experiments for testing the theory, provided 
one takes into account inevitable losses present in the 
experiment and include them into the model [ll|, [22| . 

Here we study how the field state after the click de- 
pends on detector's parameters. It is worth considering 
how SPD's actually operate: despite technical and struc- 
tural differences associated with every kind of detector, 
the photodetection process is based on the same princi- 
ple: the 'sensor' initially set in a 'ground state' interacts 
with the field and is likely to absorb a photon, doing a 
transition to the 'excited state'. After some time the sen- 
sor decays back to the ground state, emitting a photoelec- 
tron that triggers the 'amplification mechanism' (AM) 
of the SPD (e.g., by an avalanche process), producing 
a pulse of macroscopic electric current or voltage, that 
originates a registered click of the detector, representing 
one count. Besides, in real photodetectors there is a phe- 
nomenon called dark counts: photoelectrons originated 
due to the intrinsic excitations inside the AM, and not 
due to the absorption of one photon from the field. The 
influence of dark counts on the results of various experi- 
ments was considered in [H, [U, [25|, [26|, [27} , and different 
schemes of calculating dark count probability and single- 
photon quantum efficiency were used in [28l . l29l . |30| . 

In the CPM, all the results concerning the photode- 
tection process are described by means of a single entity 
characterizing the photodetector — the Quantum Jump 
Superoperator (QJS) J, that represents the back-action 
of the detector on the field upon a single photodetection. 
Immediately after the photodetection an initial field state 
described by the statistical operator p changes abruptly 
to p' = Jp/Tr[Jp], and the probability for registering 
one count during the time interval [t, t + At) is Tr[Jp]Ai, 
where At is the time resolution of the detector. It is sup- 
posed that At is small compared to other characteristic 
time scales and that QJS is time- independent, in an ideal 
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case being given by 

jp^ l0 d P d\ (i) 

where O is a lowering operator responsible for subtracting 
one photon from the field and 7 is roughly the counting 
rate [l8| with dimension (time -1 ). 

Srinivas and Davies proposed ad hoc O — a, where 
a is the bosonic lowering operator. In spite of having 
some inconsistencies as noted by the authors themselves, 
this QJS has been widely used since then. Recently, an- 

a 1/2 

other QJS defined with O = E- = (n +1) 'a (where 
n = a)a) was proposed, also ad hoc, in [3ll . |32| - we 
named it 'E-modeP for simplicity. The differences be- 
tween these two QJS's were studied in [l8j], showing that 
the inconsistencies of the SD-modcl are indeed elimi- 
nated. 

In [33| we have proposed a microscopic model (some 
other models were considered in {3~il [3q ) for the field- 
detector interaction, where we showed that both QJS's 
proposed ad hoc are particular cases that can be derived 
from a general time-dependent 'transition' 1 superopera- 
tor. However, there we considered a simplified model 
of the detector at 'zero temperature', by assuming that 
there were no intrinsic excitations inside the photodetec- 
tor and that the sensor was at exact resonance with the 
field mode. Here we relax these conditions, being less 
stringent, we take into account the effects implied by a 
'non-zero temperature' detector possessing intrinsic exci- 
tations and allow a detuning between the field and sensor 
frequencies. Moreover, we attribute numerical values for 
the model parameters in order to reproduce experimental 
data both qualitatively and quantitatively. We show that 
the dark counts appear naturally in our model, and com- 
paring the predicted counting rates and Signal-to-Noise 
ratio with experimental data, we get a good qualitative 
and quantitative agreement. Furthermore, we demon- 
strate that the actual ex pre ssion for the QJS should be 
incremented as (cf. Ref. [231 ]) 

jp = l0 d P d^ + lD D P b\ (2) 

where 7 D is roughly the dark counts rate and the oper- 
ator D describes the back-action of the detector on the 
field due to dark counts. Finally, we point out that by 
simply modifying the field frequency we are able to en- 
gineer the QJS, thus obtaining either the SD-model or 
E-model in specific regimes. 

The paper is organized as follows. In Sec. II we model 
the sensor of the photodetector as a 2-level system, ac- 
cording to the well known Jayncs-Cummings model; tak- 
ing into account the effects of the sensor-AM coupling, 
we obtain an explicit form of the transition superoper- 
ator, from which we derive a general expression for the 
QJS. In Sec. Ill we compare our results for counting rates 
with experimental data and obtain specific expressions of 
the QJS for different field wavelengths. The section IV 
contains the summary and conclusions. 



II. MODELLING THE PHOTODETECTOR 

We assume the SPD as being constituted of two parts: 
the sensor and the AM. The sensor is modelled as a two- 
level quantum object with resonant frequency loq, inter- 
acting with the mono-modal EM field with frequency u 
(for multi-modal field one should just consider a frequen- 
cies distribution). It has a ground and an excited states 
\g) and \e) (before and after the photoabsorptions), so we 
describe it by the usual Jaynes-Cummings Hamiltonian 
[3(| (for its applicability see, e.g. [13, HU and references 
therein) 

H = -Wooo +ujn + gaa+ + g*a^a- , (3) 

where g is the sensor-field coupling constant, and the 
sensor operators are cto = |e)(e| — \g)(g\, o+ = \e)(g\ 
and a- = |<?)(e|. We assume g being real, since only its 
absolute value enters the final expressions. 

Upon absorbing a photon the sensor initially in the 
ground state jumps to the excited state and some time 
later it decays back emitting a photoelectron into the 
AM. The AM is a complex macroscopic structure that 
depends on the type of SPD; it amplifies the photoelec- 
tron and originates some observable macroscopic effect, 
giving rise to the clicks of the detector. In order to de- 
scribe general features of the AM independent of the type 
of SPD, we model it as a thermal reservoir with a mean 
intrinsic excitations number n. Thus, the whole system 
field-SPD is described by the effective 'standard master 
equation' 39] 

P T = ~ [Ho, Pt\ - 7" {°-°+Pt - Zv+PtV- + Pt< j - <t +) 

i 

- 7 (n + I) ((T + (T_/3 T — 2(7_/9 T (T + + p T (T + (T_) , (4) 

where 7 is the sensor-AM coupling constant (we do not 
consider the field damping due to cavity losses). 

According to CPM the trace of QJS gives the proba- 
bility density p(t) for the photodetection, i.e. emission 
of a photoelectron at time t, given that at time t = the 
detector-field system was in the state 

Po = \9)(9\®P, (5) 

where the field state is p. Microscopically this means 
that in the time interval (0, t) the sensor has undergone 
a transition \g) — > |e) due to the absorption of one pho- 
ton, and p(t)At is the probability for decaying back to 
the ground state during the time interval [t, t+At), emit- 
ting a photoelectron that will lately originate one click. 
Following the quantum trajectories approach (39j |. p(£) 
is calculated from the master-equation ([3]) by identifying 
the sensor decay superoperator 

Rp = 2j(n+i)a-p o-+ (6) 

describing the instantaneous |e) — > \g) decay and the 
consequent emission of a photoelectron. The no-decay 
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supcropcrator UtPo = Pu(t) describes the non-unitary 
evolution of the field-SPD system from t = to t > 
without emission of photoelectrons; Pu(t) is the solution 
to the master equation ^ without the decay term ©: 



where for I > 



S,(t)p = f &!■■•[ dtihpel 
Jo Jo 



(16) 



j t Pu = -i( H ePu - Pu H l) + 2-/na +Pu a_, (7) 
where the effective non-Hcrmitian Hamiltonian is 
He = - — ~ ero+um+gao'++gcfl(T- — »7(ra+ — ). (8) 

Thus the probability density for the observation of a pho- 
tocount, or a click, at time t is equal to 



p(t) = Ttf-d RUtPo 



where RU t p represents an evolution of the field-SPD 
system from initial state p at time t = to the time t 
without any decay of the sensor and an instantaneous de- 
cay at the time t. Moreover, taking the trace only of the 
detector variables, one obtains the expression describing 
the action of the detector on the field upon the click - a 
predecessor of QJS that we call transition superoperator 



E(t)p = Trc RU t p 



(9) 



from which the probability density for a count is 
p(t)=Tr[H(t)^ 

Thus, for obtaining the transition superoperator one 
should first determine the no-decay superoperator Ut and 
then evaluate Eq. I© using the initial state §5§. In order 
to solve Eq. (J7J) we do the transformation 



Pu — XtpjjX}, 

where 

X t = exp(-iH e t), 
and obtain a simple equation for ~p v 
d 



dt 



p LJ = 2jncr + p u a-, 



(10) 

(11) 

(12) 
(13) 



~ Pu (t) = Po + 2 7 n / dt'Z + (t')pu(-t')Z-(t')- (14) 
Jo 

Now one iterates Eq. Q14[) and obtains a power expansion 
in terms of n, substitutes the result into Eq. (| 1 0[) and 
then evaluates Eq. ((9]), finally obtaining 



(7+ = X-tO+Xt, <J- = er+, 
whose formal solution is 



9i(t,t lr -- ,U) = (e\X t a + (h)a+(t2) ■ ■ -a+(ti)\g). (17) 

For I = the integrals and the terms a + should be 
dropped out: 

E (t)p = 6 p0l 6 (t) = {e\X t \g). (18) 
After some algebraic manipulations Hl[ we get for 

CD 

X t = exp[-7i(n + l/2)-iami] (19) 

x {xn+i(t)\e)(e\+Xn(-t)\g)(g\ 

-ie- luJt/2 S fl+1 (t)aa + - ie iujt ' 2 S n {t)a ] a-}, 

where 

C A (t) = cos (rftBfi/b) , S A (t) = sin (TtBfi/b) /B A , (20) 



Xn=e-^ t/2 [C A (t) - i5S A (t)] 



B;, 



Vn + 6 2 , 



(21) 



(22) 



8 = (q-ib)/2, q = (u -u)/g, b = >y/g. (23) 

As will be shown later, in realistic cases we need only 
the first three terms of the expansion (|15|) . whose con- 
stituents are found to be 

9 Q = -ie-^+V2)-Mfi+i/2)t 5fl+i ^ a (24) 

h = e-^+W-^x^t-tM-h) (25) 
9 2 = -ie-^+^-^-^Xn+iit-h) (26) 

XSfifa - t 2 ) Xn-li-h)^ . 

Substituting these expressions into Eqs. (fTS)) - l|18[). 
the transition supcropcrator turns out to be time- 
dependent, contrary to the definition of a QJS. We pro- 
ceed as in (33|, defining the QJS as the time average of 
the transition superoperator over the time T (to be deter- 
mined later) during which the photoelectron is emitted 
with high probability: 



(27) 



Considering the weak coupling, lu,loq 3> 7, \g\, under 
which both the Jayncs Cummings Hamiltonian ([3|) and 
the master equation (j4|) are valid, and expressing the field 
density operator in Fock basis as 



E(t)p = 2bg (1 + n) £ (2 7 n)< %{t)p, (15) 



P= Yl Pmn\ m ){n\ 
m,n—0 



(28) 
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after the averaging in ([27)) the off-diagonal elements of Jp 
vanish due to rapid oscillations of the terms exp(±iu)t). 
This means that the photodetection destroys the coher- 
ence of the density matrix. This can be understood from 
the point of view of information theory: information flows 
from the field-sensor system to the AM, so decoherence 
is active; moreover, since counting photons informs only 
about diagonal elements, which arc proportional to the 
number of photons, non-diagonal elements can be com- 
pletely ignored. Therefore, from now on we shall treat 
only the diagonal elements of Jp in (|2"7| . Applying the 
superoperators 5j on the density matrix as in Q15p . after 
evaluating Eq. (|27[) one is left with 



(29) 



Jp = Pnn nj( B) \n-l)(n-l\ 



n=0 



+4 D) \n)(n\ + (n + l)4 E) \n + l)(n+ 1| + • • ■ 

After some lengthy however straightforward calcula- 
tions one obtains the following expression for the first 
term in 

bg (1 +n)F n - G n 



j(B) 



\B„ 



(30) 



— (l + 2n - exp [-r(l + 2n)] 



F„ 



(1 + 2n)2 + ^ 
( [(l + 2TT)cos(r0 n )-0„sin(r^„)]|, (31) 

r 1 - exp [-2t(1 + 2n)] 

2 + Mi + m) ' lf e " = ±(1 + 2n) 



G, 



IT, 



2Re(B n ) 



otherwise, 



(32) 

b , f„ s H!^>. <3 3 , 

Since expression (|30[) is too involved to be interpreted an- 
alytically, we shall treat it numerically, as well as ongoing 
terms. In order to evaluate them we expand the time- 
dependent functions (C„ and S n , Eq. ([20]) ) in terms of 
exponentials and integrate the resulting expressions, ob- 
taining for the second term in (|29|) (here for complex B n 
defined by Eq. fl52J)) 



where 



D) _ 9 nb2 i 1 + n 



— W 3 W k V- 



j,fe=i yj y k i=i 



i 



«7 



(34) 



a x = i(l + m) + ((Jj - u* k ) /b, 

a 2 = a.\ + (% - /&, 

W 1( 2) = (l-S/B n+1 )(l-(+)6/B n ), 

W m = W m (B n+1 -> -B n+1 ) , 

lOi = w 2 = -w 3 = -w 4 = B n+1 , 

Vi = ~V4 = -{B n +i + B n ), 

V2 = —ya = B n — B n+1 . 



In the same manner one can obtain exactly all the further 
terms, but we shall not write the resulting expressions 
here. 

The QJS (|29[) contains an infinite number of terms, so 
every time the detector emits a click, the field state p 
reduces to an incoherent mixture (due to the decoher- 
ence process described above) of different states, each 
one with its respective probability. Let us examine these 
terms more closely. The first term, with coefficient J„ , 
takes out a photon from the field and modifies the rel- 
ative weight of the state components, so it represents a 
click preceded by a photoabsorption - we call this event 
a 'bright count'. The second term, dependent on 
and proportional to h (quite small as will be shown be- 
low), docs not subtract photons from the field but only 
modifies the relative weight of the state components - 
it represents the dark count, when the detector emits a 
click due to the amplification of its intrinsic excitations. 
All further terms in Eq. (|29[) are proportional to n l , 
I > 2; they describe emissions of several photons into the 
field upon a click, so we call the first of these term J„ 
the 'emission term'. We would like to stress again that 
all these processes happen simultaneously with different 
probabilities upon a click of the detector, so the post- 
measurement state of the field is a classical mixture of 
these events. We also note that there are many different 
phenomena that give rise to dark counts [j| [l2|, HI] ; our 
model takes into account only those of them that cause 
the transition of the sensor from the ground to the excited 
state, and since the sensor state depends on the sensor- 
field interaction, the dark counts modify indirectly the 
relative weight between the field components, thus de- 
pending on n. The other physical phenomena that do 
not cause this transition should not in principle modify 
the relation between the field components and should be 
expressed by a constant term equal to the corresponding 
dark counting rate. We shall not consider these effects 
since we assume that all the intrinsic excitations occur 
within the sensor. 



III. EXPERIMENTS AND QJS ENGINEERING 

Now we shall compare the predictions of our model 
with the available experimental data. Experimentally, 
the dependences of both bright and dark counting rates 
are set as functions of the wavelength of the light and the 
'bias parameter' (BP) of the detector. By BP we englobc 
such quantities as bias voltage, bias current or any other 
physical parameter the experimentalist adjusts in order 
to achieve simultaneously the best Signal-to-Noise ratio 
S and the highest bright counting rate. S is the ratio 
between the bright counting rate TZb and the dark one 
TZd, S = IZb/TZ-d, and it has the following useful prop- 
erty: as one increases the value of BP, the bright count- 
ing rates increases while the S remains the same until 
the breakdown value of BP, after which the iS starts to 
fall rapidly as function of BP. So most detectors usually 
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FIG. 1: Signal-to-Noise ratio as function of b at resonance 
(A = 500 nm) and far from it (A — 1500 nm) in the inset. 
The estimated breakdown value is bs — 380. 



FIG. 2: Bright counting rate as function as wavelength of 
the field for different values of b, showing that TZb increases 
proportionally to b. The resonant wavelength is Ao = 500 nm. 



operate near the breakdown BP in order to achieve the 
optimal performance. Experimentally TZb is deter- 
mined by directing laser pulses containing single photons 
at a given repetition rate on the detector and calculating 
the rate of counts, so in our model it is described by the 

( B) 

term J{ ; analogously, TZb is calculated in the absence 
of any input signal, so it is given by Jq D ^ . 

To do the comparison with experimental results first 
we need to set the values of our model free parameters: 
uiq, g, t and n; for the sake of better comparison we shall 
express the frequencies ujq and u> in terms of respective 
wavelengths Ao and A. Thus we are left with only two 
experimental variables: A and 6, where b plays the role 
of BP. Meanwhile, our general model can't take into ac- 
count the BP x b dependence for every kind of detector; 
nevertheless, one can argue that BP and b must be pro- 
portional to each other, since for zero BP one should 
also have 6 = 0, since in this case the detector would 
be turned off. Fortunately, we do not need to know the 
exact dependence of BP on b provided we determine the 
breakdown value bs corresponding to the breakdown BP 
at the resonance and take it as a measure of b. So bs 
will be our last fixed parameter, even though dependent 
on the free parameters, needed to compare the model 
predictions with experiment. 

After several numerical simulations we have chosen the 
following values for our model free parameters that re- 
produce qualitatively the common experimental behav- 
ior [3|, 0, H, @ and lie within the applicability region of 
the Jaynes-Cummings Hamiltonian ([3]) and the master 
equation (|4]): A = 500 nm, g = 10 11 Hz, t = 5 x 10 5 and 
n = 10 -11 , so bs — 380, as shown in figure 1. Notice that 
1) the values of 5 ~ 10 5 at the resonance (A = 500 nm, 
q = 0) and 10 2 'far away' from it (A = 1.5 /im, q 3> 1) 



are quite realistic, and 2) the chosen mean number of 
intrinsic excitations n is much larger than the number 
of thermal photons as calculated from Planck's distribu- 
tion for room temperature (rip ~ 10 -30 ), meaning that 
the contribution from defects within the detector is much 
stronger than the one from thermal photons (remember 
that we are considering a specific part of dark counts, as 
described in section II) . Moreover, we verified that below 
bs both TZb and TZb depend approximately linearly on 
&, which is in agreement with our qualitative arguments. 

In figure 2 we have a plot of the TZb for two differ- 
ent values of b as function of the field wavelength, which 
shows a good agreement, qualitatively and quantitatively, 
with experimental results and illustrates the fact that TZb 
increases proportionally to b. We could confirm numeri- 
cally that TZd does not depend on the field wavelength as 
expected, because the dark counts are detector's internal 
events. 

So our model agrees qualitatively in all aspects with 
the experimental data. Now we turn to our main goal: 
how the QJS depends on the experimental detector pa- 
rameters? First, we check that for the chosen parameters 
the emission terms and further terms in Eq. (|29p ) 

are at least 10 orders of magnitude smaller than the dark 
counts term and even more for bright counts term, which 
confirms that detectors indeed do not emit photons into 
the field. Intuitively, the photon emissions by the detec- 
tor would be possible only at temperatures much higher 
than room temperature through black body radiation, 
which is not the case in experiments. 

Thus, in practice one is dealing only with bright and 
dark counts terms that act on the field simultaneously 
every time a count is registered, so the QJS has the fol- 
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FIG. 3: Normalized bright counts term as function of n in 
di-log scale at breakdown bs for different field wavelengths. 
At resonance (A = 500 nm) we have j3 ~ 1/2 and far away 
from resonance (A = 1000 nm) j3 ~ 0. 



FIG. 4: Normalized dark counts term as function of n at 
breakdown bs for different field wavelengths and the same 
graph in di-log scale in the inset. 



lowing (diagonal) form in Fock basis 

Jp = diag (Jb + Jd)p 



(35) 



In figure 3 we show the dependence of the normalized 

( B} ( B) 

bright counts term J„ / J\ on n in di-log scale (for 
better visualization we joined the points). We note that 
near and far away from resonance one has nearly poly- 
nomial dependence (linear in di-log scale) 



J {1 



TZb n 



-20 



(36) 



with [3 ~ 1/2 at resonance and (3 ~ far away from it. 
Thus, in these cases one can write the operator depen- 
dence of bright counts term as 



JbP = 'R'B (n+1) 13 a pa) (h+1) 



(37) 



thus recovering the E-model with (3 ~ 1/2 at the reso- 
nance and SD-model with (3 ~ far away from it (A = 1 
/im). This is an important result: by just modifying the 
wavelength of the field one can engineer the QJS, obtain- 
ing one of the ad hoc proposals or another. 

Now we turn to the normalized dark counts term 
JrF^ I ' Jq D \ shown in figure 4 in linear scale and in di- 
log scale in the inset. First, we see that out of reso- 
nance (A = 1 pm) J„ is almost independent on n, so 



TZd = const. At resonance 

(D) 



we can set in this case 

(A = 500 nm), we see that for n — one gets J^' —TZd 

and for n > we have J^>' ~ ^ ' • n~ 2/3 , where 
~ 1/2 and d is a number less than 1. This means that 
at resonance the dark counts are suppressed by the pres- 
ence of light, being generated predominantly by vacuum 



fluctuations. This can be understood by the following ar- 
gument: the dark counts occur when the detector, in the 
ground state, is intrinsically excited by its excitations; 
but at resonance, the rate of excitations by field photons 
is much greater than the one by the intrinsic excitations, 
so the dark counts 'have no time' to appear and therefore 
become suppressed. Thus, the operator form of the dark 
counts term in Eq. (l29l) is 



J D p = lZ D [AopA + dAn p pn A] , 



(38) 



where Ao = |0)(0|, A = 1 — Ao and at resonance we have 
E-model with (3 ~ 1/2 and d < 1. Far away from the 
resonance we recover the SD-model with (3 ~ and d = 1. 
Thus by operating near breakdown BP and by varying 
the wavelength of the field one can engineer the QJS and 
predict its bright and dark counts terms behavior; the 
only inconvenience for obtaining SD-model is that the 
S ratio is smaller than for E-model, since in this case 
one should operate far away from the resonance. For the 
sake of completeness we could also add to the right-hand 
side of Eq. (|3"51) a constant term (IZd'p) for the dark 
counts that do not cause the transition \g) — > \e) inside 
the sensor, however the present model does not embraces 
such phenomena. 



IV. SUMMARY AND CONCLUSIONS 

We presented a microscopic model for a realistic pho- 
todetector in which we modelled it as a 2-level quan- 
tum sensor plus a macroscopic amplification mechanism. 
Using the quantum trajectories approach we deduced 
a general QJS describing the back-action of the detec- 
tor on the field upon a photocount and showed that it 
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can be represented formally as an infinite sum of terms. 
In that sum we have identified the terms corresponding 
to the bright counts (real photoabsorptions), the dark 
counts and emission events, each one occurring with its 
respective probability. Adjusting the free parameters 
of the model to fit experimental data, we showed that 
the emission terms can be disregarded in realistic sit- 
uations since their contribution becomes insignificant, 
so the QJS consists effectively only of bright and dark 
counts terms. Moreover, we reproduced qualitatively and 
quantitatively the experimental behavior of the counting 
rates and the Signal-to-Noise ratio, showing the break- 
down phenomenon. Finally, we showed that with the 
detector operating near its breakdown bias one can en- 
gineer the QJS by modifying the wavelength of the field. 



In particular, one recovers the QJS's proposed previously 
ad hoc: at resonance one gets the E- model, and far away 
from it the SD-model is identified. Last but not least, the 
contribution of the dark counts to the QJS was derived 
within the context of a photocount model. 
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